AD- A 1 33  656  A  MOLECULAR  ORBITAL  STUDY  OF  ATMOSPHERICALLY  IMPORTANT  1// 
SPECIES! U )  COLLEGE  OF  THE  HOLY  CROSS  WORCESTER  MA  DEPT 
OF  CHEMISTRY  C  A  DEAKYNE  JUN  83  AFOSR - TR -83-0803 
UNCLASSIFIED  AF0SR-82-0198  F/G  7/4  NL 


r 


▼ 


* 


AFOSR-TR-  83-  0  803 


m 


au: 


A  MOLECULAR  ORBITAL  STUDY  OF  ATMOSPHERICALLY  IMPORTANT  SPECIES 


Type  of  Report  and  Period  Covered:  Final  Report,  June  -  3  /  1983. 

Author:  Carol  A.  Deakvne 

Department  of  Chemistry 
College  of  the  Holy  Cross 
Worcester,  MA.  01610 


Controlling  Office:  United  States  Air  Force 

Air  Force  Office  of  Scientific  Research 
Building  410 

Bolling  AFB ,  D.C.  20332 


Security  Class:  Unclassified 


Distribution  Statement:  Approved  for  public  release,  distribution  unlimited 

Date  of  Submission:  August  18,  1983 


tim-P-our 


Accession  For  , 

NTTS  GHAftI  J 

DTIC  TAB  □ 

Unannounced  □ 

Justification - 


By - 

Distribution/ _  _ 

Availability  Codes 
.Avail  and/or 
Dlat  I  Special 


DTIC 


SELECTE 
OCT  171963 


D 


8  «->  Jl  o  i  i 


020 


AIR  70RCE  OFFICE  OF  SCIENTIFIC  Resrartw  MW* , 
NOTICE  OF  TTOlJfSMITTAI.  TO  -TIC  ^SE^°H 
This  technical  i*r *nn-+  v.n..  U  . 

s™r  r ii: ^  “ 


wiairiDut ion  is  unlimited,  - 

RESEARCH  OBJECTIVES  MATTHEW  J .  JCERP'ER 

Cbi.f,  TechMoal  Information  Division 

Recently  balloon-borne  mass  spectrometers  have  been  utilized  to  determine  the 


ionic  composition  of  the  earth's  ionosphere  and  stratosphere .  These  studies 

indicate  that  both  the  positive  and  negative  ions  are  clustered  species.  Consequent¬ 
ly,  the  properties  of  some  of  these  cluster  ions,  such  as  their  structures,  ener- 

8  —  11 

getics,  and  charge  distributions,  have  been  under  extensive  investigation. 

Information  on  these  properties  is  needed  for  the  following  reasons.  First,  data 

on  the  structures  of  these  complexes  will  1)  lead  to  a  more  complete  understanding 

of  the  forces  between  ions  and  neutral  molecules,  in  particular  the  role  of  hydrogen 

bonding,  which  is  believed  to  be  sigm f leant , 1 °  ’  ' ^  ’ 1  1  in  the  solvation  of  positive 

and  negative  ions  and  2)  provide  information  on  the  number  of  solvent  molecules 

which  can  be  accomodated  bv  an  ion  and  on  its  inner  a"d  outer  solvation  shells. 

This  data  is  useful  since  the  kinetic  properties,  electron  affinities,  mobilities, 

and  products  of  the  ionic  recombination  reactions  of  these  clusters  are  believed  to 

1  4 

He  dependent  on  their  sizes.  Second,  thermodynamic  data  on  clustering  reactions 
is  required  in  the  theoretical  development  of  aerosol  formation10  and  in  the  theo¬ 
retical  development  of  the  chemical  mechanisms  involving  negative  ions  which  cause 


radio  and  radar  blackout,  i.e.  photodetachment  and  photodissociation. 


Third,  the 


calculated  charge  distributions  will  indicate  whether  the  positive  or  negative 
charge  is  relatively  evenly  distributed  or  whether  it  is  localized  as  is  assumed  in 
some  electrostatic  models  used  to  compute  bond  energies. 

In  particular,  cluster  ions  of  NO^  with  H^O  have  received  considerable  attention 
recently.10’11'1^  18  Wu  and  Tiernan1  Vave  generated  ions  of  this  type  from  reaction 
mixtures  containing  various  combinations  of  NO,  NO^,  H^O,  0 2>  and  O^.  One  type  of 
product  ion  they  reported  has  the  general  formula  OH'-HNO^.  Only  one  ion  with  this 
general  formula  (bond  dissociation  energy  approximately  75  kcal/mole)  was  obtained. 
These  results  agree  with  those  of  Paulson18  and  Yamdagni  and  Kebarle.17 


One  objective  of  this  work  was  to  use  ab  initio  molecular  orbital  (MO)  theory  to 


^  H(.  =. 


I 


determine  the  optimum  smir-urc  or  oA^- HNO^pand 


to  compute  its  adiabatic  electron 


affinity,  bond  di -sedation  energy  and  charge  distribution.  Thus,  a  study  of 

Nl  7 

OH  -HNO^and  a  series  of  simpler  molecules  and  ions,  which  contain  similar  types 
of  bonds  and  whose  structures  and  adiabatic  electron  affinities  are  known,  was 

carried  out  to  find  the  most  economical  basis  set  which  yields  reliable  results. 

_  \- 

The  other  molecules  and  ions  considered  were  OH,  OH  ,  NH ,  NH~,  NH^ ,  NH^\",  00H,  00H  , 
NO^»  NOp  ,  HHd,  OH~-H^O,  HNO  ,  and  HNO^'  .  Basis  sets  with  and  without  polarization 
functions,  diffuse  functions,  and  configuration  interaction  were  utilized  to  in¬ 
vestigate  their  effect  on  the  properties  of  the  above  species. 

„  „  A  second  objective  of  this  work  was  to  probe  the  accuracy  of  the  various  basis 

sets  m  calculating  energy  changes  for  reactions  involving  negative  ions.  The  set 
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of  reactions  chosen  for  this  purpose  was  the  set  observed  bv  Paulson  and  Dale  in 
their  double  mass  spectrometer  and  selected  ion  flow  tube  (SIFT)  study  of  OH  -H^O 
and  NO.J.  Three  of  the  reactions  thev  observed  were: 


OH  -H20  +  N02 - >  HN03  +  H20 


N02  +  H20  +  OH 


N02  +■  OH  +  H20 


•Vhi, 


Th«.r  results  -  also  indicated  that  once  the  HNO^t  was  formed  it  Drake  down  via 

V 

electron  detachment  or  dissociation  into  N0'2  and  OH  around  interaction  energies 
of  0.8eV.  The  dissociation  of  HNOtjj\  into  OH**  and  NO^  was  concluded  to  be  a  less 
likely  loss  process^ 


HN03  +  e 


N02  +  OH 


N02  +  OH 


X- 


These  reactions  were  chosen  since  thermodynamic  data  is  available  for  all  of  them 
and  they  involve  manv  of  the  species  considered  above. 

COMPUTATIONAL  DETAILS 

The  calculations  were  carried  out  ab  initio  usmq  the  Guassian  80  computer  pro- 
gram^°  on  a  DEC  VAX  11/780  computer.  Mol ler-Plesset  perturbation  theory^  termi¬ 
nated  at  the  second  (MP2)  and  third  (MP3,  for  OH  and  OH  only)  o'-der  was  applied  to 
compute  electron  correlation  effects.  for  OH  and  OH  corrections  to  the  Hartree 
Fock  (HF)  energy  were  computed  also  utilizmq  a  configuration  interaction  method 
including  all  double  substitutions  (CID).  The  unrestricted  MP  method  was  employed 
for  open-shell  systems  ;the  restricted  MP  method  was  employed  for  closed-shell  systems. 
The  frozen  core  approximation,  where  onlv  valence-shell  orbitals  are  correlated,  was 
used  throughout  the  calculations. 

Several  different  basis  sets  were  utilized  in  the  work.  The  4- 31G  basis  is  an 

22 

s,p  split  valence  basis  set.  6-31G*  is  similar  but  includes  also  polarization 

9  j 

functions  of  d  type  on  atoms  other  than  hydrogen.  The  6-31G*  basis  set  augmented 

24 

with  p-functions  on  the  hydrogens  is  denoted  6-31G**;  6-311G**  is  a  triple-valence 

s,p,d  basis  set.^^  When  diffuse  sp  functions  on  atoms  other  than  hydrogen  are  in- 

26 

eluded  m  the  basis  set,  the  symbol  +•  appears  in  its  designation,  i.e.,  4-31+G. 

Standard  scale  factors  and  exponents  were  employed  for  all  basis  functions  other 

than  the  diffuse  functions.  The  scale  factor  for  the  diffuse  functions  was  1.0. 

The  exponent  for  the  diffuse  functions  added  to  nitrogen  was  0.053  and  to  oxygen  was 
26 

0.068.  These  exponents  yielded  better  results  than  those  suggested  by  Dunning  and 

27 

Hay  (0.048  for  nitrogen  and  0.059  for  oxygen).  For  molecules  containing  both 
nitrogen  and  oxygen,  calculations  were  carried  out  with  and  without  diffuse  functions 
on  the  nitrogen. 

The  geometries  of  the  molecules  and  ions  were  optimized  via  the  Fletcher-Powell 
28 

method  to  minimize  the  total  energies.  Reported  bond  lengths  represent  convergence 


to  0.001  X  and  bond  angles  to  0.1°.  The  total  energies  were  then  used  to  compute 
electron  affinities  (E.A.)  and  bond  dissociation  energies  (D)  according  to  the 
following  equations. 

E.A. (B)  =  E(B)  -  E(B-)  (7) 

DUB)  =  E(B)  +  E  ( A )  -  E(AB)  (8) 

No  zero-point  corrections  are  made  for  any  of  the  calculated  energies  unless  other¬ 
wise  stated. 

STATUS  OF  RESEARCH 

OH,  OH  .  The  structures  and  energies  of  OH  and  OH  and  the  adiabatic  electron 
affinity  of  OH  obtained  at  the  various  basis  set  levels  utilized  in  this  work  are 
presented  in  Tables  I  and  II.  These  data  are  compared  with  the  experimental  and 
with  other  theoretical  results  in  the  tables.  The  structures  of  these  species  and 
of  all  the  other  species  considered  in  this  study  are  optimized  at  the  basis  set 
noted  unless  stated  otherwise.  The  following  abbreviations  are  employed  in  the 
tables:  DZP,  a  basis  set  of  double-zeta  quality  which  includes  polarization  functions; 
CEPA,  coupled  electron  pair  approach  to  the  treatment  of  electron  correlation;  PNO-CI, 
configuration  interaction  treatment  based  on  pseudo  natural  orbital  configuration  ex¬ 
pansion;  EOM,  equations-of -motion  technique  for  calculating  electron  affinities; 

CGTO,  contracted  Gaussian  type  orbitals. 

The  results  show  the  following:  1)  Including  correlation  in  the  basis  set 
lengthens  the  O-H  bond  lengths  while  including  polarization  functions  shortens  them. 
They  are  shortened  more  when  polarization  functions  on  hydrogen  are  included  than 
when  they  are  not.  2)  When  diffuse  functions  are  added  to  the  split-valence  s,p 
basis,  the  O-H  bonds  are  shortened;  however,  when  they  are  added  to  an  s,p,d  basis, 
the  O-H  bond  lengths  are  essentially  unchanged.  3)  All  of  the  split-valence  HF  optimi- 
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TABLE  I.  Theoretical 
of  OH. 

and  Experimental 

Structures,  Energies  and 

Election  Affinities 

Method 

ROH 

(X) 

Energy 
( hartrees ) 

Electron  Affinity 
(eV) 

MNDOa 

0.937 

— 

0.27 

HF/STO-3Gb 

1.014 

-74.36489 

-8.16 

HF/3-21Gb 

0.986 

-74.97023 

-2.77 

RHF/DZPC 

0.958 

-75.39980 

1.03 

PNO-CI/DZP+d 

— 

— 

1.27e 

CEPA/DZP+d 

— 

— 

1.5le 

SCF-CI/DZPf 

0 . 970 ( assumed ) 

— 

1.19 

SCF-FOM/DZPf 

0 .970 ( assumed) 

— 

1.46 

HF/6-3 lG*b 

0.959 

-75.38228 

-1.52 

MP2/6-3lG*g 

0.979 

-75.52103 

-0.22 

MP3/6-31G*g 

0.981 

-75.53300 

-0.44 

CID/6-31G*h 

0.980 

-75.53085 

-0.52 

MP2/6-311G»*h 

0.967 

-75.57290 

0.019 

HF/4-31+Gh 

0.955 

-72.77317 

0.090 

MP2/4-31+Gh 

0.986 

-75.38242 

1.61 

MP2/6-31+Gh 

0.991 

-75.46104 

1.58 

HF/6-31+G*h 

0.959 

-75.38705 

-0.26 

MP2/6-31+G*h 

0.981 

-75.52931 

1.68 

MP3/6-31+G*h 

0.982 

-75.54069 

1.14 

CID/6-31+G*h 

0.981 

-75.53633 

1.11 

MP2/6-31+G**h 

0.973 

-75.54050 

1.74 

MP2/6-31+G**(5D)h 

0.974 

-75.53862 

1.73 

Expt. 

0.9711 

-75 . 780  3 

1.83* 

aReference  29.  bReference  30.  CReference  31.  dReference  32.  SCorrected  for  zero-point 
energies.  ^Reference  33.  ^Reference  34.  bThis  work.  ^"Reference  35.  -'Reference  36. 
^Reference  37. 
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TABLE  II.  Theoretical  and  Experimental  Structures  and  Energies  of  OH  . 


Method 

ROH 

<8. 

Energy 
(hartrees ) 

MNDO3 

0.939 

— 

HF/STO-3Gb 

1.068 

-74.06502 

HF/3-2lGb 

1.029 

-74.86863 

RHF/DZP° 

0.955 

-75.36650 

PNO-CI/DZP+d 

0.956 

-75.67423 

CEPA/DZP+d 

0.961 

-75.69325 

HF/fe- 3lG*b 

0.962 

-75.32660 

MP2/6-31G*0 

0.981 

-75.51314 

MPi/6- ilG*6 

0.981 

-75.51703 

CID/6-31G*0 

0.980 

-75.51165 

MP2/6-311G**® 

0.968 

-75.57360 

HF/4-31<-G6 

0.951 

-72.77648 

MP2/4-31+G® 

0.988 

-75.44139 

MP2/6-31+Ge 

0.993 

-75.51892 

HF/6-31+G*® 

0.953 

-75.37758 

MP2/6-31+G*® 

0.978 

-75.59099 

MP3/6-31+G*6 

0.972 

-75.58250 

CID/6-31+G*6 

0.972 

-75.57701 

MP2/6-31+G**® 

0.971 

-75.60427 

MP2/6-31+G**(5D)e 

0.971 

-75.60197 

Expt . 

0.970f 

-75.847g 

aReference  29.  ^Reference  30.  CReference  31.  ^Reference  32.  eThis  work. 
f Reference  35.  ^Reference  36. 
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zations  yield  bond  lengths  that  are  too  short  compared  to  the  experimental  values. 

4)  All  of  the  MP  and  CID  values  are  too  long,  with  the  exception  of  the  case  where 
correlation,  diffuse  functions,  and  polarization  functions  on  all  atoms  are  incorp¬ 
orated  in  the  calculations,  i.e.  MP2/6-31+G**.  The  latter  values  are  within  the 
experimental  uncertainties  of  the  observed  values.  5)  The  MP2/4-31+G  bond  lengths 
are  too  long  by  approximately  0.02?.  6)  The  MP2  technique  yielded  geometries  that 

are  as  good  or  better  and  electron  affinities  that  are  better  than  those  obtained 
from  the  MP3  or  CID  techniques.  7)  There  is  not  necessarily  a  correlation  between 
more  stable  total  energy  and  more  accurate  electron  affinities.  8)  Adding  diffuse 
functions  to  an  s,p,d  split-valence  basis  set  improves  the  accuracy  of  the  calculated 
electron  affinities  approximately  as  much  as  adding  electron  correlation  does.  9) 

All  of  the  calculated  electron  affinities  are  too  low  implying  that  the  negative  ion 
is  not  represented  well  enough  and  that  zero-point  corrections  are  important  for 
these  systems.  10)  For  this  work,  only  the  electron  affinity  of  OH  computed  at  the 
MP2/6-31+G**  and  MP2/4-31+G  levels  are  reasonable.  These  values  differ  from  the  ex¬ 
perimental  value  by  O.leV  and  0.2eV,  respectively.  11)  Polarization  functions  are 
not  necessary  to  attain  reasonable  electron  affinities.  12)  The  electron  affinities 
obtained  from  this  work  are  as  good  or  better  than  other  values  obtained  from  theo¬ 
retical  calculations  reported  in  the  literature.  13)  No  improvement  in  the  results 
is  achieved  by  using  the  MP2/6-31+G  basis  rather  than  the  MP2/4-31+G  basis,  by  ex¬ 
pressing  the  polarization  functions  in  term  of  six  d-orbitals  rather  than  five,  or 
by  employing  a  triple-valence  basis  set  (MP2/6-311G** ) . 

Atomic  charges  from  population  analysis  for  OH  and  OH-  and  dipole  moments  for 
OH  at  several  basis  set  levels  are  given  in  Tables  III  and  IV,  respectively.  For  OH, 
all  of  the  calculated  charge  separations  are  too  large  since  the  dipole  moments  are 
too  large.  This  is  expected  to  be  true  for  OH-  and  all  of  the  other  species  studied 
as  well.  The  agreement  between  the  theoretical  and  experimental  dipole  moments  is 
improved  by  adding  polarization  functions  to  the  basis. 
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TABLE  III.  Atomic  Charges 

and  Dipole  Moment 

for  OH 

Method 

Charge  on 

Charge  on 

Dipole  Moment 

0 

H 

(debyes  ) 

HF/6- 3lG*a 

-0.444 

0.444 

1.88 

MP2/6-31G*b 

-0.442 

0.442 

1.90 

MP2/6- 31 1G**C 

-0.252 

0.252 

1.85 

MP2/4-31+GC 

-0.408 

0.408 

2.23 

MP2/6-31+GC 

-0.416 

0.416 

2.27 

HF/6-31+G*° 

-0.473 

0.473 

1.96 

MP2/6-31+G*C 

-0.471 

0.471 

2.01 

MP2/6-31+G**C 

-0.364 

0.364 

1.96 

d 

Expt . 

— 

— 

1.66 

aReference  30.  ^Reference  34.  CThis  work.  ^Reference  39. 


TABLE  IV.  Atomic  Charges  for  OH 


Method 

Charge  on  0 

Charge  on  H 

HF/6-31G*3 

-1.206 

0.206 

MP2/6-31G*b 

-1.200 

0.200 

MP2/6-311G**C 

-1.063 

0.063 

MP2/4-31+GC 

-1.269 

0.269 

MP2/6-31+GC 

-1.282 

0.282 

HF/6-31+G*C 

-1.351 

0.351 

MP2/6-31+G*C 

-1.343 

0.343 

MP2/6-31+G**C 

-1.241 

0.241 
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Based  on  these  results,  it  was  decided  to  consider  only  two  basis  sets,  where 
feasible,  for  the  species  remaining  to  be  investigated,  i.e.  MP2/4-31+G  and  MP2/6-31+G** . 
The  latter  basis  includes  five  d-orbitals. 

NH,  NH-.  Structures  and  energies  of  NH  and  NH-  and  adiabatic  electron  affinities 
of  NH  are  given  in  Tables  V  and  VI.  The  charge  distributions  obtained  from  the  MP2/ 
4-31+G  and  MP2/6-31+G  **  (in  parentheses) calculations  are  given  below. 


(-0.243) 

N 

-0.234 


(0.243) 

H 

0.234 


(-1.136)  (0.136) 

N  —  H 
-1.123  0.123 


The  dipole  moment  of  NH  computed  at  the  MP2/6-31+G**  level  is  1.82D.  Several  obser¬ 
vations  can  be  made  upon  examining  the  data  in  Tables  V  and  VI  and  in  the  above  fig¬ 
ures.  1)  The  MP2/4-31+G  bond  length  is  about  O.OlX  too  long  for  NH  and  0.02?  too  long 
for  NH  compared  to  their  experimental  values.  2)  The  MP2/6-31+G**  bond  length  is 
within  experimental  error  for  NH  and  about  0.01&  too  long  for  NH- .  3)  All  of  the 

calculated  electron  affinities  are  too  low,  with  the  best  value  from  this  work  too 
low  by  0.5eV.  The  MP2/6-31+G**  and  MP2/4-31+G  values  vary  by  about  O.leV.  4)  Since 
it  is  unlikely  that  zero-point  corrections  to  the  energies  will  improve  the  theoretical 
results  by  O.SeV,  a  more  accurate  basis  set  is  required  to  properly  describe  the 
electron  delocalization  in  NH-.  5)  The  above  results  imply  that  the  charge  separations 
from  population  analysis  are  too  large.  6)  The  polarization  functions  have  very  little 
effect  on  the  charge  distributions  of  these  systems. 

NHj ,  NH^-.  Bond  lengths,  bond  angles,  and  total  energies  of  NH^  and  NH2  and 
adiabatic  electron  affinities  of  NH^  are  tablulated  in  Tables  VII  and  VIII.  MP2/ 
4-31+G  and  MP2/6-31+G**  (in  parentheses)  charge  distributions  are  shown  below. 


(0.536)  -0.599 


H  H 

(0.268)  0.299 


(-1.340)  „  -1.392 

H  H 


(0.170) 


0.196 
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TABLE  V.  Theoretical  and  Experimental  Structures,  Energies  and  Electron  Affinities 
of  NH. 


Method 

rnh 

d) 

Energy 
( hartrees ) 

Electron  Affinities 
(eV) 

MNDOa 

0.993 

-0 . 7 8 

HF,/DZP+b 

— 

— 

-1.43C 

b 

PNO-CI/DZP+ 

— 

— 

-0. 25C 

b 

CEPA/PZP+ 

— 

— 

0 . 0 1 c 

d 

MP2/4- 31+G 

1.046 

-54.942S6 

-0.220 

d 

MP2/6-31+G** 

1.035 

-55.06980 

-0.151 

Expt . 

1 .036e 

— 

0.  38* 

3 

Reference  29. 

b 

Ret erence 

32. 

°Corrected 

for  zero-point  energies. 

d 

This  work. 

e 

Reterence  35. 

^Reference 

40. 

TABLE  VI.  Theoretical  and 

Experimental  Structures  and 

Eneraies  of  NH  . 

Method 

rnh 

(8) 

Energy 
( hartrees ) 

o 

a 

1.000 

— 

HF/DZP+b 

1.022 

-54.92219 

PNO-CI/DZP+b 

1.033 

-55.12025 

CEPA/DZP+b 

1.039 

-55.13703 

MP2/4-31+GC 

1.058 

-54.93448 

MP2/6-31+G**C 

1.044 

-55.06426 

d 

Expt. 

1.037 

— 

aRef erence  29.  bReference 

c  d 

32.  This  work.  Reference 

41. 
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TABLE  VII.  Theoretical  and  Experimental  Structures,  Energies,  and  Electron 
Affinities  of  NH^. 


Method 

rnh 

(8) 

e 

H  NH 
(°) 

Energy 
( hartrees ) 

Electron  Affinity 
(eV) 

CNDOa 

1.016 

105.4 

— 

0.61 

b 

MNDO 

1.002 

104.4 

— 

-0 . 44 

CGTO/DZP+C 

MP2/4-31+Gd 

1.016 

105.4 

— 

0.61 

1 .029 

108.8 

-55.56300 

0.387 

MP2/6-31+G**d 

1.025 

e 

104.0 

-55 .71336 

0 . 464 

f 

Expt . 

1.024 

103.3 

0 . 779“ 

Reference  42. 

^Reference  45. 

b  c 

Reference  29.  Reference 

43.  This 

e 

work.  Reference  44. 

TABLE  VIII.  Theoretical 

and  Experimental 

Structures 

and  Energies  of  . 

Me  thod 

rnh 

\H 

Energy 

( A  ) 

<°> 

(hartrees ) 

MNDO3 

1.013 

100.9 

— 

CGTO/DZP+b 

1.016 

103.9 

— 

MP2/4-31+GC 

1.037 

107.7 

-55.57728 

MP2/6-31+G**C 

1.029 

103.7 

-55.73041 

Expt.  (solid)d 

1.03 

104 

— 

aReference  29.  Reference  44.  CThis  work.  ^Reference  46. 
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TABLE  VII.  Theoret  1  cal  and  Kxrfrimonr.il  Structures,  Energies,  and  Electron 
Affinities  of  NH  ,  . 


Method 

rnh 

($, 

G 

H  NH 
(°) 

Energv 
( hartrees ) 

Electron  Affinitv 
(eV) 

CNDO3 

1 . 016 

105 . 4 

— 

0.61 

h 

MNDO 

1.002 

104.4 

— 

-0 . 44 

cgto/dzp+c 

1.01b 

105.4 

— 

0.61 

MP2/4-  11+g'J 

1 . 029 

108.3 

-55.56300 

0  .  387 

MP2  /6-31-eG**d 

1 . 025 

104.0 

-55.71336 

0.464 

Kxpt . 

1.0246 

103.3® 

— 

f 

0  .  77<T 

a 

Reference  42. 

^Reference  45. 

b  c 

Reference  29.  Reference 

43.  dThls 

e 

work.  Reference  44. 

TABLE  VIII.  Theoretical 

and  Experimental 

Structures 

and  Energies  of  NH^  . 

Method 

rnh 

©H 

NH 

Energy 

(8) 

t°) 

(hartrees ) 

MNDO  3 

1.013 

100.9 

— 

CGTO/DZP+b 

1.016 

103.9 

— 

MP2/4-31+GC 

1.037 

107.7 

-55.57728 

MP2/6-31+G**C 

1.029 

103.7 

-55.73041 

Expt.  (solid)d 

1.03 

104 

— 

aReference  29.  ^Reference  44.  CThis  work.  ^Reference  46. 
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The  MP2/6- 31-t-G**  dipole  moment  of  is  2.21D. 

The  following  conclusions  can  be  drawn  from  the  tables  and  figures.  1)  The 

N-H  bond  lengths  are  reproduced  verv  well  at  both  basis  set  levels  for  NH0  and  NH^  . 

2)  These  basis  sets  vield  bond  lengths  closer  to  the  experimental  values  than  the 

other  theoretical  results  reported.  3)  The  agreement  between  the  experimental  and 

calculated  bond  angles  is  also  verv  good  at  the  MP2/6-31+G**  level  for  both  species. 

However,  the  MP2/4-3WG  bond  angles  are  too  large.  4)  The  best  electron  affinity 

obtained  from  this  work  is  approximately  0.3eV  lower  than  the  experimental  value  and 

O.lSeV  lower  than  the  best  theoretical  value.  The  two  electron  affinities  computed 

fr^m  this  research  vary  bv  about  O.leV.  5)  These  eaicualted  electron  affinities  are 

0.2oV  closer  to  the  experimental  result  than  thev  are  for  \'H. 

H^O.  Optimized  structures  and  total  energies  for  H^O  are  presented  in  Table  IX 

Charge  distributions  from  population  analysis  are  given  in  the  figures  below  for  the 

MP2/4-31+G  and  MP2/6-31+G**  (in  parentheses)  basis  sets.  The  experimental  dipole 

39 

moment  for  H^O  is  1.95D  and  the  MP2/6-31+G**  value  is  2.28D. 

(-0.741)  -0.923 

O 

/  \ 

H  H 

(0.370)  0.461 

TABLE  IX.  Theoretical  and  Experimental  Structures  and  Energies  of  H^O. 


Method 

roh 

(X) 

®HOH 

(°) 

Energy 
(har trees ) 

CNDO/2a 

1.04 

104.5 

-19.892 

RHF/DZPb 

0.947 

105.2 

-76.0410 

HF/4-31++GC 

0.956 

106.6 

-75.94946 

MP2/4-31+Gd 

0.969 

111.7 

-76.04841 

d 

MP2/6-31+G** 

0.964 

105.5 

-76.23042 

Expt . 6 

0.957 

104.5 

aReference  47.  b Reference  31.  CReference  48.  The  4-31G  basis  set  was  augmented 

with  two  uncontracted  diffuses  functions  on  each  hydrogen  and  nine  sets  of  uncon- 

d  © 

tracted  diffuse  s  and  p  functions  on  oxygen.  This  work.  Reference  49. 
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The  results  show:  1)  The  MP2/4-31+G  and  MP2/6-3 1  <-G*»  optimized  0-H  bond 
lengths  are  about  0.01  $  too  long  and  the  optimized  H-O-H  bond  angle  is  about  7° 
too  large  for  the  former  basis  set  computation.  2)  The  calculated  charge  sepa¬ 
rations  are  too  large. 

OOH,  OOH  .  Geometries  and  energies  of  OOH  and  OOH  and  the  adiabatic  electron 
affinity  of  OOH  are  shown  in  Tables  X  and  XI.  Atomic  charges  obtained  from  the 
MP2/4-31+G  and  MP2/6-31+G**  (in  parentheses)  calculations  are  displayed  below.  The 
computed  dipole  moment  ( MP2/6- 31+G* * >  of  OOH  is  2.15D. 

(-0.7281  (-0.588)  (0.316) 

0  -  O  -  H 

-0.628  -0.733  0.361 


(-0.0297)  (-0. 363)  (0. 392) 

O  -  O  -  H 

0.00200  -0.471  0.469 


From  the  data  in  the  tables  and  figures  it  is  evident  that:  1)  The  MP2/6-31+G** 
and  experimental  bond  lengths  and  bond  angles  for  OOH  are  verv  close,  whereas  the 
MP2/4-31+G  values  are  too  large  bv  0.05$  (0-0),  0.01$  (O-H)  and  2°.  2)  For  OOH  , 

the  MP2/4-31+G  calculated  0-0  distance  varies  from  the  MP2/6-31+G**  distance  bv 
0.1$  and  the  O-H  distances  vary  by  0.015$.  The  O-O-H  bond  angles  are  similar.  3) 

It  is  necessary  to  include  polarization  functions  in  the  basis  set  for  a  proper 
description  of  the  0-0  bond.  4)  This  is  the  first  case  where  the  electron  affinity 
is  too  high  (by  0.3eV)  when  d-functions  are  not  included  in  the  basis  set  and  too 
low  (by  O.leV)  when  they  are.  Incorporating  polarization  functions  stabilizes  the 
neutral  OOH  molecule  more  than  the  OOH  ion  which  decreases  the  magnitude  of  the 
theoretically  determined  electron  affinity.  For  all  of  the  other  species  considered 
above,  the  anion  was  stabilized  more  than  the  neutral  system.  5)  The  negative 
charge  is  delocalized  more  at  the  MP2/4-31+G  level  of  calculation  but  for  both 
levels  most  of  the  additional  negative  charge  in  OOH  is  on  the  terminal  oxygen 
(0.6  -  0. 7e) . 
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TABLE  X.  Theoretical  and  Experimental  Structures,  Energies,  and  Electron 
Atfmities  of  OOH. 


Method 

Roo 

(X) 

roh 

(8) 

^OOH 

<°> 

Energy 
( hartrees ) 

Electron 

Affinities 

(eV) 

INDOa 

1  . 185 

1 .04  7 

110.6 

— 

— 

CI/DZb 

1.458 

0.973 

104.6 

-150.2448 

— 

MP2/4- 31+GC 

1  .  385 

0.985 

106.9 

-150.19329 

1.49 

MF2/6-  31.G**C 

1 . 329 

0.977 

104.7 

-150.51427 

1.06 

d 

d 

d 

O 

Expt . 

1.335 

0.977 

104.1 

1  .  19' 

a Reference  30. 

bReference  51. 

c 

This  work. 

"^Reference  44. 

e 

Reference 

52. 

TABLE  XI.  Theoretical 

and  Experimental 

Structures 

and  Energies 

of  OOH~ 

Method 

Roo 

roh 

^  OOH 

Energy 

(S) 

<8> 

(°) 

(hartrees ) 

HF/STO- 3Ga 

1.449 

1.021 

105.4 

-147.93655 

MP2/4-31+Gb 

1.610 

0.980 

96.6 

-150.24786 

MP2/6-31+G**b 

1.516 

0.966 

97.3 

-150.55336 

aReference  30.  bThis  work. 
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OH  -HjO.  Table  XII  gives  optimized  structures  and  total  energies  of  OH-H^O  . 
One  point  of  interest  with  respect  to  the  geometry  of  this  ion  is  whether  it  is 

symmetric  or  not.  Earlier  studies  have  shown  that  calculations  carried  out  with  a 

475354  5455 

limited  basis  set  * '  '  and/or  only  partial  geometry  optimization  '  are  not 

sufficiently  accurate  to  settle  this  point.  Later  work^  undertaken  bv  performing 

a  global  geometry  optimization  with  a  correlated  wavefunction  indicates  that  this 

anion  is  symmetric.  However,  no  diffuse  functions  were  included  in  the  basis  set 

of  any  of  these  prior  calculations,  and  diffuse  functions  are  known  to  be  important 

26  27 

for  u  prorer  description  of  anionic  species.  ' The  preliminary  MP2/4-31+G 

fully  optimized  results  imply  that  cis,  trans,  and  gauche  OH  -H.,0  are  asymmetric 

with  an  0-0  bond  length  of  about  2.6$  and  with  the  central  hydrogen  moving  about 

0.08A  along  the  0-0  bond  axis.  However,  the  structural  data  obtained  for  00H  and 

OOH  suggest  that  these  bond  lengths  are  too  long.  Since  it  is  widely  accepted 

now  that  a  hydrogen  bond  will  be  symmetric  when  the  distance  between  the  heavv  atoms 
o  58 

is  less  than  2.50a  and  the  corrected  0-0  bond  distance  will  be  close  to  this 

value,  an  even  more  accurate  basis  set  (i.e.,  MP2/6- 31 +G*» )  is  required  to  resolve 

this  question.  In  addition,  other  conformations  must  be  considered  since  Rohlfing 

et.al.56  find  that  the  dihedral  angle  defined  by  the  external  hydrogen  is  110°. 

At  the  MP2/4-31-t-G  level  of  calculation,  the  trans  and  gauche  systems  are  of 

essentially  equal  energy,  and  they  are  more  stable  than  the  cis  system  by  1.70  kcal. 

All  three  systems  have  slightly  non-linear  hydrogen  bonds  with  the  central  hydrogen 

5.3°,  6.1°,  and  1.9°  above  the  0-0  bond  axis  in  the  gauche,  trans,  and  cis  forms, 

respectively.  The  bonds  between  the  oxygens  and  terminal  hydrogens  have  remained 

essentially  constant  in  length  (to  within  0.01$)  compared  to  their  values  in  OH 

and  H20  and  the  H-O-H  angle  has  decreased  by  about  3°.  These  results  are  consistent 

with  Rohlfing  et  al’s.  work^  and  with  results  from  other  studies  of  asymmetric 
58 


hydrogen  bonds. 
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The  adiabatic  electron  affinity  of  OH-H^O  can  be  computed  from  the  MP2/4-il+G 

total  energies  of  OH  -H^O,  OH,  and  H^O  bv  assuming  that  the  bond  dissociation  energy 

(D)  of  OH-H^O  is  negligible.  The  number  obtained  is  2.86eV  which  is  in  good  agree- 

19 

ment  with  the  experimental  value  of  2.91eV.  The  latter  is  derived  from  E.A.  (OH) 
and  D  (OH  -H.,0)  assuming  D  (OH.H_>C)  is  negligible. 

A  comparison  of  the  charge  distributions  calculated  at  the  MP2/4-31-I-G  level 
for  OH  ,  H?0  and  trans  OH  -H^O  (given  below)  shows  that  there  is  a  transfer  of  only 
0.0500e  from  the  OH  to  the  H^O  upon  formation  of  the  hydrogen  bond.  The  electron 
density  is  transferred  to  the  oxygen  and  terminal  hydrogen  of  the  water.  However, 
most  of  the  gain  in  electron  density  on  these  two  atoms  comes  from  the  central  hydro¬ 
gen  which  loses  0.129e  when  the  hydrogen  bond  is  formed.  Similar  charge  redis¬ 
tributions  were  found  for  the  cis  and  gauche  structures  and  in  other  studies  of 

58 

hvdrogen  bonded  systems. 


H 

0.389 


O  — 

■1.026 


0.590 

H 


xv 


329 


-1.282 


NO^#  NC>2  .  Structures  and  energies  of  NO^  and  NO^  and  electron  affinities  of 
NC^  are  tabulated  in  Tables  XIII  and  XIV.  From  these  tables  it  is  observed  that: 

1)  For  NO^  the  MP2/6-31+G**  and  experimental  NO  bond  lengths  differ  by  0.025$  and 
ONO  angles  differ  by  0.4°.  The  distance  and  angle  computed  at  the  MP2/4-31+G  level 
are  too  large  by  0.6$  and  0.4°,  respectively.  2)  For  NO^  the  MP2/6-31+G**  and  gas 
phase  experimental  NO  bond  lengths  vary  by  0.13$  and  the  ONO  angles  vary  by  3.5°. 
The  distance  calculated  at  the  MP2/4-31+G  level  is  too  large  by  0.18$  and  the 
angle  is  too  small  by  4°.  3)  In  all  other  cases  where  a  comparison  can  be  made, 

the  MP2/6-31+G**  and  experimental  bond  lengths  differ  by  no  more  than  0.04$.  This 
suggests  that  the  experimentally  determined  gas  phase  NO  distance  may  be  too  short 


TABLE  XII.  Theoretical  and  Experimental  Structures  and  Energies  of  OH  ‘H^O. 


Method 

a 

r, 

a 

r2 

a 

r  3 

a 

r4 

<Pa 

a 

Energy 

<8> 

(8) 

(8) 

(8) 

<°> 

<°> 

(°» 

<°) 

( hartrees  ) 

CNDO/2(sym)b,'_ 

1.04 

1.17 

1.17 

1.04 

106 

106 

180.0 

180.0 

-  38.928 

HF/ST0-3G(asvm)C,d 

1.007 

1.010 

1  .  199 

1  .  163 

100.9 

101 . 1 

180.0 

180.0 

-149. 14322 

HF/4-3lG(asym)C,d 

0.955 

1.106 

1  .  348 

0 . 9b2 

108.8 

108.0 

180.0 

180.0 

-151 .20343 

HF,’4-31G(  svm  )6  ’ ' 

0.958 

1.212 

1.212 

0.958 

109.2 

109.2 

180.0 

150.0 

-151 .20325 

HF/4-31G(asym)e’^ 

0.955 

1.097 

1  .  361 

0.962 

108.1 

110.0 

176.1 

135.0 

-151 .20358 

c  d 

HF/6- 31G* ( asvm )  ’ 

0.948 

1.095 

1.  384 

0.951 

102.6 

101.7 

180.0 

180.0 

-151 . 39246 

HF/6-31G*(svm)e,t 

0.950 

1.207 

1.209 

0.950 

102.2 

10  2.2 

180.0 

120.0 

-151 . 39176 

0  f 

HF/6-31G*,,asvm)  ’ 

0.948 

1 . 040 

1.488 

0.953 

100.8 

101.8 

173.2 

115.0 

-151 . 39340 

HF/6-3lG**(asym)d’^ 

0.949 

1.095 

1.384 

0.951 

102.6 

101.7 

180.0 

180.0 

-151.41148 

HF/6-31G**( svm )e ’ ^ 

0.94b 

1.20  1 

1.206 

0.946 

102.7 

102.7 

180.0 

12  0. 0 

-151.41105 

HF/6-31G**( asvm)^’ t 

0.944 

1.044 

1.457 

0.949 

101.4 

102.5 

1  74  .  5 

115.0 

-151 .41228 

Cl  'DZ ( sym ) b  ’ 1 

0.958 

1.206 

1 . 206 

0 . 958 

104.5 

104.5 

180.0 

180.0 

-151. 72689 

CI/DPZ ( asvm )h ’  1 

0.957 

1.092 

1  .  373 

0.95  7 

104.5 

104.5 

180  .  1 

180.0 

-151.91518 

MP2/6-3lG**(svm)e’^ 

0.965 

1.223 

1  _  f  2  J 

0 . 9b4 

99. 1 

99. 1 

180.0 

no.  o 

-151.82012 

MP2/4-31+G(asvm)k’1 

0.971 

1.047 

1.56  t 

0.979 

108.2 

128.3 

168.3 

180  .  ) 

-151.5 3584 

MP2/4- 31 +G( asvm ) k ’ 1 

0.971 

1.048 

1.558 

0.980 

108.0 

123.9 

171.1 

60.0 

-151 .53585 

MP2/4-31+G( asvm ) k ’ 1 
Expt . m  (solid) 

0.971 

1.040 

1.145 

1.586 

1.145 

0.980 

108.3 

126.0 

179.8 

0.0 

-151.53313 

3  be 

See  figure  below  for  definitions.  Reference  47.  The  hydrogen  bond  was  assumed  to  be 

d  e 

linear  and  trans.  All  other  geometrical  parameters  were  optimized.  Reference  53.  All 
geometrical  parameters  were  optimized.  ^Reference  56.  ^The  energy  was  calculated  using 
the  HF/6-31G*  optimum  structure.  hThe  hydrogen  bond  was  assumed  to  be  linear.  Only  r ^ 
and  r_  were  optimized.  The  other  parameters  are  taken  from  experimental  values  for  H.,0. 

i  i  i  k  1 

Reference  54.  JReference  55.  This  work.  The  dihedral  angles  were  fixed  at  the  values 
given  in  the  table.  All  other  parameters  were  optimized.  Reference  57.  The  hvdrogen 
bond  was  assumed  to  be  symmetric. 


0^  is  the  dihedral  angle  derived  from  the  external  hydrogens. 
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bv  about  0.10$.  4)  Baird  and  Tavlor  find  that  adding  polarization  functions  to 

onlv  the  nitrogen  shortens  the  N-O  bonds.  This  result  implies  that  optimized  geo¬ 
metries  obtained  using  the  MP2/4-31+-G  basis  should  be  improved  by  adding  polari¬ 
zation  functions  to  the  nitrogen.  It  is  also  possible  that  this  addition  will  im¬ 
prove  the  theoretical  electron  affinities.  (It  is  being  tested  now.)  5)  Including 
diffuse  functions  on  N  in  the  basis  set  has  little  effect  on  the  optimized  geometries 
but  does  improve  the  electron  affinities.  Thus,  it  is  recommended  that  structures 
be  determined  with  a  basis  set  which  does  not  incorporate  diffuse  functions  on 
nitrogen  and  then  that  these  structures  be  used  to  compute  a  new  total  energy  with 
a  basis  which  does  incorporate  these  functions.  6)  The  best  MP2/6-31+G**  value  for 
the  electron  affinity  is  too  low  by  0 . 25-0 . 35eV ;  the  best  MP2/4-31+G  value  is  too 
low  bv  0.05-0. lOeV.  Again  the  electron  affinity  calculated  without  polarization 
functions  is  larger  ( bv  0.2eV)  than  that  calculated  with  polarization  functions  in 
the  basis  set.  Although  a  more  complete  analysis  is  required,  it  is  possible  that 
there  will  be  a  compensatorv  effect  between  neglecting  polarization  functions  and 
zero-point  energies  for  systems  containing  more  than  one  non-hvdrogen  atom.  7)  The 
best  electron  affinity  obtained  from  this  work  is  as  good  or  better  than  the  numbers 
reported  in  the  literature. 

Table  XV  lists  charge  distributions  of  NO?  and  NO,,-  from  population  analysis 
and  dipole  moments  of  NO^  It  is  clear  from  these  data  that:  1)  The  calculated 
dipole  moments  are  too  large.  2)  Incorporating  diffuse  functions  on  N  in  the  basis 
set  decreases  the  charge  separation  in  the  N-0  bonds  and  slightly  improves  the 
dipole  moments.  3)  Although  more  of  the  extra  negative  charge  on  NO?  is  distributed 
to  the  oxycpns  than  the  nitrogen,  the  percentage  of  electron  density  delocalized  on 
the  nitrogen  increases  to  30%  when  polarization  and  diffuse  functions  on  N  are  in¬ 
corporated  in  the  basis  set.  4)  This  is  the  only  set  of  molecules  studied  thus  far 
where  there  is  a  significantly  greater  charge  separation  in  the  N-0  bonds  and  a 
slightly  poorer  computed  dipole  moment  when  polarization  functions  are  included  in 
the  basis  than  when  they  are  not  included. 


TABLE  XIII.  Theoretical 
Af  f ini  ties 

and  Experimental  STructures,  Energies, 
of  NOn 

and  Electron 

Method 

rno 

^ONO 

Energy 

Electron  Affinity 

(&> 

(°) 

( har trees ) 

( eV ) 

MNDO3 

1 .  174 

133.0 

— 

2.55 

HF/4-31Gb 

1.19 

136.7 

-203.69429 

1.51 

N*  b 

HF/4- 31G 

1  .  16 

135 

-203.77780 

1 .11 

SCF-EOM/DZ+C 

1.21 

134 

-203.9609 

2.25d 

P  .  f 

MP2/4- 31+G 

1.250 

133.7 

-204.10858 

2.14 

MP2/4- 3l+Ge ,q 

1.250 

133.7 

-204.11125 

2.21 

MP2 /6- 3 1+G**e ’ f 

1.217 

133.7 

-204.56441 

1.93 

MP2/6-31<-G**e,g 

1.217 

133.7 

-204.56579 

2.02 

Expt . 

1 . 193h 

1 34.  lh 

— 

2 . 275 1 

Expt . 

— 

— 

— 

2.  -61 

■3  b 

Reference  29.  Reference  59.  The  N*  designates  that  polarization  functions  were 

c  d 

added  to  the  nitrogen  only.  Reference  60.  Corrected  for  zero-point  energies. 

e  f  g 

This  work.  Diffuse  functions  were  added  to  the  oxygen  only.  'Diffuse  functions 

were  added  to  the  oxygens  and  nitrogen.  The  energy  was  calculated  at  the  optimized 

geometry  obtained  without  diffuse  functions  on  the  nitrogen.  ^Reference  61. 

Reference  62.  ^Reference  63. 

TABLE  XIV.  Theoretical  and  Experimental  Structures  and  Energies  of  NO^  . 


Method 

rno 

^ONO 

Energy 

(8> 

(°) 

( hartrees ) 

MNDOa 

1.215 

116.5 

HF/4-3lGb 

1.256 

117.0 

-203.74991 

N*  b 

HF/4-31G 

1.23 

116 

-203.81844 

MP2/4-31+GC,d 

1.328 

115.6 

-204.18704 

MP2/4-31+GC,e 

1.328 

115.6 

-204.19229 

MP2/4-31+GC,f 

1.327 

115.8 

-204.19229 

MP2/6-31+G**C,d 

1.280 

115.7 

-204.63531 

MP2/6-31+G**C,e 

1.280 

115.7 

-204.64003 

MP2/6-31+G**C,f 

1.278 

116.0 

-204.64005 

Expt.  (gas)g 

1.15 

119.5 

Expt.  (crystal )b 

1.236 

115.4 

See  footnotes  on  the  following  page. 


aReterence  29.  ^Reference  59.  CThis  work.  ^Diffuse  functions  were  added  to  the 
0 

oxvqens  only.  Diffuse  functions  were  added  to  the  oxygens  and  nitrogen.  The 
energy  was  calculated  at  the  optimized  geometry  obtained  without  diffuse  functions 
on  the  nitrogen.  ^Diffuse  functions  were  added  to  the  oxvgens  and  nitrogen. 
^Reference  62.  ^Reference  64. 


HNChj,  HNO^  .  Table  XVI  gives  optimized  geometries  and  total  energies  of  HNO, 
and  planar  HNO,  and  adiabatic  electron  affinities  of  HNC  .  .  Pvramidal  structures 
of  UNO,  are  currently  under  investigation,  and  the  preliminary  data  indicate  that 
pvramidal  HNO ^  is  the  more  stable  form.  The  MP2/4-31+G  charge  distributions  for 
HNO,  and  HNO^  are  presented  below.  The  calculated  dipole  moment  of  HNO ^  is  2.99D; 
the  experimental  value  is  2.190. 


-0.411  O  0.497 

I 

N  0.522 

/  \ 

0  o 

-0.269  -0.339 


-0.595  0  0.439 

I 


N  0.404 


/  \ 

O  0 

-0.1 98  -0.650 


The  following  conclusions  can  be  drawn  from  these  results.  1)  The  optimized 
bond  distances  between  the  nitrogen  and  terminal  oxygens  differ  from  the  experimental 
numbers  by  0.04-0.05$.  The  other  calculated  N-0  distance  is  too  long  by  about  0.10$. 
The  bond  angles  are  within  2°  of  each  other.  2)  Based  on  earlier  results,  the  com¬ 
puted  O-H  bond  is  likely  to  be  about  0.015$  too  long.  3)  Some  of  the  variance  in 
calculated  and  observed  geometries  is  due  to  experimental  uncertainty.  4)  At  present 
the  theoretically  obtained  electron  affinity  is  0.4eV  too  low.  This  error  will  be 
improved  when  the  most  stable  form  of  HNO^-  is  determined  and  when  diffuse  functions 
on  the  nitrogen  are  incorporated  in  the  basis  set.  It  is  expected  that  the  error 
in  the  computed  electron  affinity  will  then  be  comparable  to  that  found  for  NO., . 


TABLE  XV 

Atomic  Charges  and  Dipole 

Moments  for  NO 

2  and  N02  . 

Molecule 

Method 

Charge 

Charge 

Dipole  Moment 

on  N 

on  0 

(debyes ) 

NO 

MP2/4-31+Ca 

0.548 

-0.274 

0.687 

C 

b 

-0.0945 

0.672 

MP2/4-31+G 

0.189 

MP2/6-31+G**3 

0.659 

-0.329 

0.748 

b 

MP2/6-31+G** 

0.386 

-0.193 

0.733 

„  c 

Expt . 

— 

— 

0.316 

NO," 

MP2/4-31+Ga 

0.325 

-0.662 

c 

b 

MP2/4-31+G 

-0.0303 

-0.485 

MP2/4-3UGd 

-0.0276 

-0.486 

MP2/6-31+G**a 

0.403 

-0.701 

MP2/6-31+G**b 

0.0843 

-0.542 

MP2/6- 31+G** 

0.0838 

-0.542 

aDif fuse 

functions  were  added  to  the 

oxygens  only. 

^Diffuse  functions  were  added 

to  the  oxygens  and  nitrogen.  The  energy  was  calculated  at  the 

optimized  geometry 

obtained 

without  diffuse  functions  on 

the  nitrogen. 

c 

Reference 

39.  ^Diffuse 

functions  were  added  to  the  oxygens  and  nitrogen. 


5)  The  charge  distribution  and  structure  of  HNO^  enables  one  to  answer  the  question 
as  to  whether  it  should  more  properly  be  considered  OH  -NO^.  The  atomic  charges 
show  that  HNO^-  is  the  more  correct  designation  since  the  OH  moiety  has  a  total 
charge  of  -0.15e.  Furthermore,  the  fact  that  the  bond  between  the  nitrogen  and 
the  oxvgen  of  the  OH  group  is  only  0.025$  longer  in  HNO^  than  in  HNO^  argues  that 
the  anion  is  a  covalent  compound  rather  than  a  complex.  6)  Every  atom  in  HNO^  gains 
electron  density  when  it  is  negatively  charged,  but  most  of  the  extra  electron  den¬ 
sity  is  distributed  to  the  terminal  oxvgens  ( 0 . 3e  each).  7)  As  usual  the  theoretical 
determined  dipole  moment  is  too  large. 

H_NO.  .  Earlier  work  performed  in  this  laboratorv  led  to  the  conclusion  that 

2  4 

OH  -HNO^  cannot  be  considered  a  solvated  point  charge  and  is  more  correctly  desig¬ 
nated  as  H^NO^  .  Several  possible  conformers  of  this  species  have  been  investigated 
at  several  basis  set  levels  (Table  XVII).  The  MP2/4-31+G  calculation  remains  to  be 
done.  A  complete  geometry  optimization  was  carried  out,  with  the  exception  of  the 
0-H  bond  lengths  in  conformer  II  (HF/6-31G*  calculation)  and  conformer  III  (HF/4-31G 
calculation).  The  earlier  study  revealed  that  the  effect  on  the  total  energy  of  the 
system  of  optimizing  the  0-H  bond  lengths  is  negligible  compared  to  the  effects  of 
the  other  parameters.  All  three  conformers  have  a  distorted  tetrahedral  geometry. 
Both  the  HF/4-31G  and  HF/6-31G*  bond  lengths  are  most  likely  too  short,  and  the 
bond  angles  are  most  likely  too  large. 

From  the  total  energies  given  in  Table  XVIII,  it  is  evident  that  conformer  I 
is  only  slightly  more  stable  than  the  other  two  conformers.  The  bond  dissociation 
energy,  D(H2NC>4  ),  is  calculated  according  to  the  following  equation  and  is  about 
0 . 87t:V  too  low  compared  to  the  experimental  number. 

D(H.NO  ")  =  E(HNO  )  +  E(OH_)  -  E(H,NO  ~) 

2  4  3  2  4 

The  discrepancy  in  these  values  may  be  due  in  part  to  experimental  uncertainty,  since 
Ddl^NO^-)  wa8  never  measured  directly  but  was  determined  from  a  thermodynamic  cycle. 
In  addition,  if  any  conformation  not  yet  considered  is  more  stable  than  I,  the  mag- 
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TABLE  XVI.  Theoretical  and  Experimental  Structures,  Energies,  and  Electron 
Affinities  of  HNO^  and  HNO ^  . 


Molecule 

Method 

a 

ri 

(X> 

a 

r2 

(X) 

a 

r  3 

<x> 

a 

r4 

(X) 

a 

© 

<°> 

a 

<s> 

o 

f  ) 

a 

(°» 

Energy 
( har trees ) 

Electron 

Af f in i tv 

( eV ) 

HNO. 

3 

HF/4-31Gb 

1 . 202 

1 .  202 

1.421 

0.960 

115.2 

122.4 

100.0 

-278.98916 

_ - 

( planar ) 

HF/4-31GC 

1.219 

1  .  194 

1  .  373 

0.961 

128.9 

116.2 

107.7 

-278.99211 

— 

HF/6- 31G*' 

1.188 

1.172 

1.334 

0.955 

129.2 

116.1 

105.3 

-279.44426 

— 

MP2/4-3UG° 

1.260 

1.250 

1.50  6 

0.989 

131.1 

115.9 

104.0 

-279.54618 

0.15 

d 

Expt . 

1.21 

1.21 

1.41 

0 . 96 

130 

114 

102 

0 . 5  7e 

( assumed ) 

0.56f 

hnc>3~ 

( planar ) 

MP2/4-31+G 

1.348 

1.311 

1.530 

0.985 

132.0 

112.9 

99.0 

-279.551-4 

— 

3  b 

See  figure  below  for  definitions.  Reference  65.  r  and  r7  were  constrained  to  be  identical. 

C  *”  (j 

This  work.  Diffuse  functions  were  added  to  the  oxvaens  only.  The  structural  data  is  from 

&  f 

reference  66.  and  r ^  were  assumed  to  be  identical.  Reference  67.  Reference  19. 


nitude  of  the  theoretically  computed  DfH^NO^  )  would  increase. 

Adding  polarization  functions  to  the  basis  set  worsens  the  agreement  between 

the  theoretically  and  experimentally  obtained  bond  dissociation  energies.  Taking 

correlation  into  account  improves  the  agreement  between  them.  Thus,  the  polarization 

functions  stabilize  the  fragments  mere  than  H  NO.  while  the  reverse  is  true  for 

2  4 

the  correlation  effects. 

HF  6-31G*  atomic  charges  are  displayed  for  HNO^  and  eonformer  I  of  H  NO  in 

j  2  4 

the  figure  below.  Several  observations  can  be  made  upon  examining  this  figure. 

1 )  Q.811e  has  been  transferred  from  the  OH  to  the  HNO^  group  in  the  ion.  2)  The 
negative  charge  has  been  distributed  over  all  of  the  atoms  in  the  HNO^  moietv; 
each  atom  is  more  negative  in  the  ion  than  it  is  in  HNO,. 


/ 


0  -0.474 

I 

N  0.920 

\  .  H  0.493 

0  ' 


-0.412  -0.527 


0 

-0.694 


0  -0.694 

I 

N  °>Z67  •  H  °-«0 

/  v  </  \ 

“  ^  0  -0.609 


O  -ad 

-0.609 


.60 

H  0.420 


'rtiermodynamic  data.  Table  XIX  lists  thermodynamic  data  for  several  reactions 
involving  some  of  the  species  studied.  Neither  zero-point  nor  termperature  correct¬ 
ions  have  been  made  in  the  ^ E  values,  unless  otherwise  indicated.  Where  possible 
the/^E's  have  been  compared  at  several  different  basis  set  levels.  The  results 
show  that  1)  The  A  E's  for  the  reactions  involving  HNO^  do  not  agree  well  with 
the  A  H's  for  these  reactions,  since  the  structure  of  HNO^-  has  not  been  completely 
optimized  yet.  2)  The  A  E's,  computed  from  the  MP2/4-31+G  data,  for  the  other 
reactions  are  all  within  0.2eV  of  the^H's.  3)  The  variance  between  the  &  H ' s  and 
A  E's  calculated  at  the  MP2/6-31+G**  level  is  larger,  again  suggesting  that  there 


TABLE  XVII.  Theoretical  Structure  of  HnNO^ 


b 

Parameters 

H  ,N0 

hong 

4_  ’  11 

v°4~. 

N-O 

1.343 

(1.284) 

1  .  375 

(1.307) 

1  .  342 

N-O, 

1.343 

(1.284) 

1  .  314 

(  1.264  ) 

1.3  36 

n-o3 

1.457 

(1.430) 

1.459 

(1.431  ) 

1.455 

N-O 

1.457 

(1.430) 

1.459 

(1.431) 

1.455 

0.-H- 

0.957 

(0.950) 

0.957 

( 0.949  i 

0.95  7 

3  0 

O  -H . 

0.95  7 

( 0.950 ) 

0.95  7 

(0.9491 

0.957 

4  b 

Oj  NO, 

120.5 

( 120.2  ) 

119.8 

1119.6) 

116.9 

°1N03 

109.2 

( 108.8  ) 

110.9 

(110.1) 

110.2 

V°4 

107.1 

(107.6) 

110.9 

(110.1  ) 

110.2 

NO  ,H_ 
i  ) 

100.9 

(99. 2  ) 

100.9 

(99.  .n 

101.1 

NO  H. 

101.0 

(99.2) 

100.9 

(99.0) 

101.1 

4  b 

-0,  0.,N0, 

i  .  i 

124.5 

(124.8) 

12  3.4 

(124.0) 

126.0 

°1°2N04 

2  34.6 

(-124.8) 

2  36.5 

(236.0) 

234.0 

0.  N'O.H. 

-12.5 

( -16.0  ) 

0.4 

(0.0) 

59.7 

1  3  : 

0 ,  NO , H  . 

214.6 

(212.7) 

-0.4 

(0.0  ' 

300.  3 

aBond  lengths  are  in  3,  bend  angles  are  m  degrees.  The  values  in  parentheses  are 

the  HF/6-31G*  optimized  values;  the  others  are  the  HF/4-31G  values.  See  the  figur 

c 

The  geometry  of  conformer  III  was  not  optimized  using  the 


below  for  definitions 
HF/6-31G*  basis  set  since  it  is  the  least  stable  of  the  three 


N 


.  -1 

C 

t  M 


TABLE  XVIII. 


Total  Energy  and  Bond  Dissociation  Energy  of  H_,NO, 

2  4 


a.u.  and  E(OH")  m  -75.51286  a.u.,  i.e.  their  MP2/6-31G*//HF/6-31G*  values. 
17a. 


Reference 


Svstem 

Method 

E 

tot 

( hartrees ) 

D( H.NO  ” ) 

2  4 

(eV) 

h2n°4-. 

I 

HF/4-31G 

-354.30674 

2.31 

HF/6-31G* 

-354.82468 

1.47 

MP2/6-31G*//HF/6-31G*a 

-355.76110 

2.55b 

h2N04"> 

II 

HF/4-31G 

-354.30377 

HF/6-31G* 

-354.82241 

h2N04-, 

III 

HF/4-31G 

-354.30067 

c 

Expt . 

3.5 

aThe  total 

energy  was  calculated  at  the 

MP2/6-31G* 

level  using  the  optimized 

geometry 

obtained 

at 

the  HF/6-31G*  level.  bThis 

number  was 

computed  using  ECHNO^)  = 

-280.15453 

J 


6- 


TADLE  VIX.  Thermodynamic  Data 


React  ion 

Method 

&E(eV/mole ) 

&H(eV/mole  ) 

(for  E) 

1  . 

OH~ '  HnO  ♦  N02  HNC>3”  +  H70 

MP2/4-31+Ga 

1.21 

0.23 (estimated ) b 

o 

OH*-H  0  +  N07 — *  NO,”  +  H00  + 

OH  MP2/4-31<-Ga 

0.72 

d 

0.63 

MP2/4-  )1+GC 

0.65 

MP2/6-31+G**'’ 

0.97 

MP2/6-31+G**C 

0.88 

3  . 

oh"-h?o  — >  oh”  +  h7o 

MP2/4-31+G 

1.25(1.11)* 

1 ,08f 

MP2/6- 31 +G** 

1. 17(1. 03)e 

HF/6- 3 1G** / /HF/6- 

3lG*g  1 . 50 (  1 . 37  )e 

HF/4-31+Gh 

1.19 

4. 

hno3"  no2  +  oh" 

MP2/4- 31+Ga 

0.049 

0 . 87 ( estimated ) b 

5. 

HN03"  — ^  N02”  +  OH 

MP2/4-  31-hGa 

-0.48 

0 . 43 ( est imated )b 

6. 

OH"  +  N02  — *?  N02"  +  OH 

MP2/4-31+Ga 

-0.5  3 

-0 . 44 ( estimated )b 

MP2/4-31+GC 

-0.59 

MP2/6-31+G**9 

-0.20 

MP2/6-31+G**C 

-0.29 

aDiffuse  functions  have  been  added  to  the  oxygens  only.  J.F.  Paulson,  personal  communi¬ 
cation.  These  values  were  not  measured  directlv  but  were  calculated  from  a  thermodynamic 

c  d 

cycle.  Diffuse  functions  have  been  added  to  the  oxygens  and  nitrogens.  Reference  19. 

©  f 

The  values  in  parentheses  have  been  corrected  for  zero-point  energies.  Reference  17b. 

^Reference  53.  b  R.R. Squires ,  personal  communication  to  J.F.  Paulson. 


is  some  cancellation  of  errors  when  zero-point  enerqies  and  polarization  functions 
are  neglected.  4)  The  results  are  particularly  qood  for  D(OH  ‘HO),  which  is  en¬ 
couraging  since  bond  dissociation  energies  are  often  poorly  represented  at  lower 

,  .  ,  10,26,58 

levels  of  theory. 

SUMMARY 

1)  At  the  MP2/4-31+G  level,  all  calculated  bond  lengths  are  too  long  and  all 
calculated  bond  angles  are  too  large  (with  the  exception  of  NO^  and  NO?  )  compared 
to  their  experimental  values  and  their  values  obtained  with  bigger  basis  sets.  XH, 
where  X  =  N,  O,  bond  lenqths  are  too  long  by  0.01  -  0.02$,  with  the  difference  beinq 
closer  to  0.01$  for  the  neutral  systems  and  to  0.02$  for  the  negative  ions.  XX  and 
XY  bond  lengths  are  too  lonq  by  0.05  -  0.10$,  aqam  the  smaller  error  is  associated 
with  the  neutral  systems  and  the  larger  error  with  the  anions.  The  discrepancies 

in  the  bond  angles  range  from  0.5°  to  7°.  The  worst  disagreements  occur  for  the  XH0 
systems.  Poor  experimental  data  for  some  of  these  species  accounts,  at  least 
partially,  for  the  large  magnitude  of  some  of  these  error  ranges. 

2)  Except  for  NO?  and  NO^  ,  the  MP2/6-31+G**  optimized  bond  distances  and 
angles  are  within  0.007$  and  0.5°  respectively,  of  the  observed  values.  Again,  poor 
experimental  data  mav  explain  the  poor  theoretical  results  for  NO^  . 

3)  Electron  affinities  and  other  thermodynamics  data  computed  at  the  MP2/6-31+G* 
and  MP2/4-31+G  basis  set  levels  are  in  reasonable  agreement  with  the  experimental 
numbers.  For  the  electron  affinities,  the  biggest  disagreement  is  for  NH,  0.5eV. 
Apparently,  one  set  of  diffuse  functions  on  the  nitrogen  is  not  sufficient  to  obtain 
accurate  results  for  NH.  However,  as  the  size  of  the  nitrogen-containing  species 
increases  and  the  electron  density  is  more  delocalized,  the  agreement  between  the 
theoretical  and  experimental  values  improves.  Overall,  the  calculated  electron 
affinities  are  as  good  or  better  than  most  of  the  other  calculated  values  reported 

in  the  literature. 


4)  Even  though  it  appears  that  the  errors  in  the  theoretical  thermodynamic 


data  arising  from  neglecting  zero-point  energies  may  be  partially  canceled  by  the 
errors  arising  from  neglecting  polarizing  functions,  a  complete  assessment  of  the 
usefulness  of  these  calculations  in  predicting  thermodynamic  data  cannot  be  made 
until  zero-point  energy  effects  have  been  evaluated. 

5)  Adding  polarization  functions  to  the  basis  set  gives  rise  to  shorter  opti¬ 
mized  bond  lengths  and  smaller  optimized  bond  angles,  even  when  the  polarization 
functions  are  incorporated  onlv  on  the  central  atom.  Taking  correlation  effects 
into  account  produces  longer  optimized  bond  lengths  and  smaller  optimized  bond 
angles.  (See  references  30  and  34,  also.)  Including  diffuse  functions  in  the 
basis  has  little  effect  on  either  the  bond  distances  or  angles.  In  contrast,  all 
three  of  these  factors  lead  to  improve  theoretical  electron  affinities  and  reaction 
energies.  Combining  these  results  suggests  that  the  geometries  of  these  species 
should  be  calculated  utilizing  either  the  MP2/4-31G  or  MP2/4-31G*  basis  set,  with 
polarization  functions  on  the  central  atom  only.  Then  the  total  energies  should 

be  computed  at  the  MP2/4-31+G  or  preferably,  where  possible,  at  the  MP2/6-31+G** 
level  using  the  previously  determined  optimized  geometries. 

6)  The  preliminary  MP2/4-31+G  results  indicate  that  OH  -H^O  is  asymmetric  but 

additional  work  using  this  and  the  MP2/6-31+G**  basis  is  required  to  settle  this 

question.  It  is  possible  that  the  asymmetric  structure  is  a  local  minimum  on  the 

potential  energy  surface  rather  than  the  true  minimum  and  that  the  relative  energies 

of  the  symmetric  and  asymmetric  forms  of  OH  ‘H^O  may  change  when  the  more  accurate 

56 

calculations  are  carried  out. 

7)  Although  the  final  optimized  geometry  of  HNO^-  has  not  been  determined  yet, 
HNO^  ,  unlike  HNO^,  has  a  pyramidal  structure.  The  thermodynamic  data  for  reactions 
involving  HNO^  will  improve  when  its  optimum  geometry  and  total  energy  are  obtained. 

8)  H2N04  has  a  <^stortec*  tetrahedral  geometry.  Several  conformers  may  very 
similar  energies.  Further  work  necessary  to  improve  the  A.E's  for  reactions 


this  anion  undergoes. 


-29- 


9)  Where  a  comparison  can  be  made,  all  of  the  calculated  dipole  moments  are 
too  large  by  0.3  -  0.8D  compared  to  the  experimental  values.  This  indicates  that 
the  computed  charge  separations  are  too  large. 

FOLLOW-ON  RESEARCH 

Specific  recommendations  for  follow-on  research  include: 

1)  Optimizing  the  structures  of  all  the  species  considered  in  the  present 
study  utilizing  the  MP2/4-31G  and  MP2/4-31G*  basis  sets  to  asses  their  reliability 
compared  to  that  of  the  MP2/6-3L+G**  and  MP2/4-31+G  basis  sets.  If  the  former  two 
basis  sets  do  yield  reliable  results,  considerable  computer  time  could  be  saved  in 
carrying  out  further  studies  of  this  type. 

2)  Calculating  the  zero-point  energies  of  all  the  species  considered  in  the 
present  study.  This  will  be  done  using  the  new  GAUSSIAN  82  program  which  should  be 
available  shortly. 

3)  Using  isodesmic  reactions,  i.e.  reactions  where  the  number  of  bonds  of  a 

68 

given  type  are  retained,  to  compute  electron  affinities  and  other  thermodynamic 

data.  This  type  of  approach  often  reduces  inaccuracies  in  calulated  reaction 
68,69 

energies . 

4)  Continuing  the  investigation  of  the  structure  of  OH  -H^O  to  determine 
whether  the  hydrogen  bond  is  symmetric  or  asymmetric  and  to  compare  the  results  at 
various  levels  of  theory.  This  will  entail  finding  the  optimum  geometry  of  OH  -H^O 
at  the  MP2/4-31+G  and  MP2/6-31+G**  levels. 

5)  Completing  the  MP2/4-31+G  optimization  of  pyramidal  HNO^-.  Then  this  geo¬ 
metry  will  be  used,  i.e.  no  additional  optimizations  will  be  performed,  to  find  the 
effect  on  the  total  energy  of  HNO^  of  adding  polarization  and  diffuse  functions  to 
the  nitrogen. 

6)  Employing  the  most  accurate  basis  set  possible,  within  feasible  economic 

and  time  constraints,  to  determine  the  otpimal  structure  and  total  energy  of  H.NO  “ 

2  4 
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and  H^NO^.  These  data  would  be  used  then  to  obtain  bond  dissociation  energies  and 
electron  affinities  for  these  species. 

7)  Extending  this  research  to  include  other  anions,  in  particular  other 
nitrogen  oxides. 
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20.  ABSTRACT  (Continue  on  reverae  aide  If  neceeeery  and  Identify  by  block  number) 

The  theoretical  research  conducted  under  this  program  was  concerned  with 
a  study  of  the  structures,  total  energies,  and  charge  distributions  of  a  series 
of  atmospherically  important  species.  These  included  OH,  OH~,  NH,  NH~ ,  H2o, 
NH2«  NH2~'  °°h*  °°h~»  OH-.H2o,  N02,  N02",  HN03,  HN03~,  and  HjNO^".  The  calcu¬ 
lations  were  carried  out  ab  initio  primarily  at  the  MP2/6-31+G**  and  MP2/4-31+G 
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basis  set  levels.  The  total  energies  obtained  were  used  to  compute  energy 
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changes  for  several  reactions  involving  the  above  species  and  electron  affin¬ 
ities  for  OH,  NH,  NH2,  00H,  0H-H20,  N02 ,  and  HN03.  The  theoretical  bond 
lengths,  bond  angles,  electron  affinities,  and  reaction  energies  are  in 
reasonable  agreement  with  their  experimental  values,  particularly  for  the 
MP2/6-31+G**  results. 

An  analysis  was  performed  of  the  effect  of  diffuse  functions,  polari¬ 
zation  functions,  and  electron  correlation  on  the  above  molecular  and  thermo¬ 
dynamic  properties.  It  was  found  that  reasonable  geometries  are  obtained 
without  adding  diffuse  functions  to  the  basis  set  and  reasonable  energetics 
are  obtained  without  adding  polarization  functions. 
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